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ABSTRACT 

An  ad  hoc  method  to  solve  boundary  value  problems  which  are  posed  on 
infinite  intervals  is  to  reduce  the  infinite  interval  to  a  finite  but  large 
one  and  to  impose  additional  boundary  conditions  at  the  far  end.  These 
boundary  conditions  should  be  posed  in  a  way  so  that  they  express  the 
asymptotic  behaviour  of  the  actual  solution  well.  In  this  paper  a  rigorous 
theory  is  derived  which  defines  classes  of  appropriate  additional  boundary 
conditions.  Appropriate  is  to  be  understood  in  the  sense  that  the  solutions 
of  the  approximate  problems  converge  to  the  actual  solution  of  the  'infinite* 
problem  as  the  length  of  the  finite  interval  tends  to  infinity.  Moreover 
boundary  conditions  which  produce  convergence  with  the  largest  expectable 
order  are  devised. 
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SIGNIFICANCE  AND  EXPLANATION 


A  boundary  value  problem  on  an  infinite  interval  consists  of  a  system  of 
ordinary  differential  equations,  some  boundary  conditions  at  a  finite  point 
and  a  continuity  condition  at  infinity,  for  example  it  is  required  that  the 
solution  converge  to  a  finite  limit  as  the  independent  variable  converges  to 
infinity.  This  condition  is  problematic  when  solutions  are  sought 
computationally.  Therefore  it  is  useful  to  cut  the  infinite  interval  at  a 
point  which  is  far  out  and  to  impose  some  suitable,  so  called  asymptotic, 
boundary  conditions  at  that  far  end  and  to  solve  the  resulting  two-point 
boundary  value  problem  which  is  now  posed  on  a  finite  but  large  interval  by 
any  appropriate  code.  Difficulties  occur  in  finding  appropriate  asymptotic 
boundary  conditions  in  the  sense  that  the  solutions  of  the  approximating  two- 
point  boundary  value  problems  converge  to  the  'infinite*  solution  as  the 
length  of  the  interval  converges  to  infinity.  This  paper  devises  suitable 


asymptotic  boundary  conditions  which  produce  a  fast  -  in  most  practical 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


A  THEORY  FOR  THE  APPROXIMATION  OF  SOLUTIONS  OF  BOUNDARY  VALUE  PROBLEMS 

ON  INFINITE  INTERVALS 
Peter  A.  Markovich 


1.  Introduction 

Boundary  value  problems  on  infinite  intervals,  which  are  posed  in  the  following  way 

(1.1)  y*  =  t“f(t,y),  1  <  t  <  »,  a  e  N  u  {0}  , 

(1.2)  y  e  C([1,«l)  :  <==>  y  e  C([1,“H  and  lim  y(t)  exists  , 

t+aj 

(1.3)  b( y( 1 ) )  *  0 

where  f  :  Rn+1  +  Rn  are  often  solved  numerically  by  restricting  the  infinite  interval  to 
a  finite  but  large  one  and  by  imposing  additional  suitable  boundary  conditions  at  the  right 
end.  The  resulting  two-point  boundary  value  problem  has  the  following  forms 

(1.4)  x^  -  t“f(t,xT),  1  <  t  <  T,  T  »  1  , 

(1.5)  b(xT(1))  =  0  . 

(1.6)  S(xt(T),T)  =  0  , 

and  can  be  solved  by  any  appropriate  code.  The  questions  this  paper  answers  are  the 
following: 

1)  What  class  of  asymptotic  boundary  conditions  S(xT(T),T)  =  0  imply  convergence  in  the 
following  sense 

d.7)  «xT  -  y * 1 1 ,T]  *  0  as  T  +  ” 

where  Izl,  :=  sup  |z(t)|  and 

ta,b]  te(a,b) 

2)  which  asymptotic  boundary  conditions  yield  a  reasonably  fast  order  of  convergence. 

It  will  be  shown  that  the  admissible  boundary  conditions  have  to  be  constructed  with 
regard  to  the  invariant  subspaces  and  eigenspaces  of  the  matrix 
d.8)  Ao(yJ  != 
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and  that  the  order  of  convergence  defends  on  the  decay  properties  of  the  solution  y  of 
(1*1) ,  (1.2)  and  (1.3)  and  on  the  largest  algebraic  multiplicity  of  the  eigenvalues  of 

(1.8)  which  have  real  part  zero.  Calling  this  number  r  we  will  show  that  for  all 
admissible  boundary  conditions  (1.6) 

(1.9)  lx  -  y!  <  const.  T( a+ 1  * r ( l nT )* ^ I S( y( T ) ,T) I ,  0  <  j  <  n 

T  ;  1  t  T  J 

holds  if  the  solution  y  of  the  problem  (1.1),  (1.2),  (1.3)  is  isolated  and  decays 
sufficiently  fast. 

Work  on  the  solvability  and  asymptotic  behavior  of  the  solutions,  of  problems  of  the 
kind  (1.1),  (1.2)  and  (1.3)  .a^  ‘'een  done  by  M.  Lentini  and  H.  B.  Keller  (1980),  F.  de  Hoog 
and  Richard  Weiss  (1980a, b)  and  P.  Markowich  ( 1980a, b). 

The  paper  is  organized  as  follows:  in  Chapter  2  linear  constant  coefficient  problems 
are  treated.  Chapter  3  is  concerned  with  linear  time  varying  problems  which  have  the 
property  that  the  matrix  describing  the  system  evaluated  at  infinity  has  distinct 
eigenvalues,  in  Chapter  4  this  assumption  ls  neglected.  Chapter  5  deals  with  nonlinear 
problems  and  a  practical  problem  from  fluid  dynamics  is  dealt  with  in  Chapter  6. 

The  novelity  of  this  paper  in  comparison  to  the  above  mentioned  ones  is  that  no  severe 
assumptions  on  f^(t,y(°°))  are  made  while  Lentini  and  Keller  (1980)  require  a  certain 
convergence  behaviour  of  this  matrix  as  t  approaches  infinity. 

The  used  techniques  are  similar  to  those  used  by  de  Hoog  and  Weiss  (1980b)  who  treated 
the  case  where  f  (“,y(“))  does  not  have  an  eigenvalue  with  a  zero  real  part. 


2.  Linear  Constant  Coefficient  Problems 


The  problem 

(2.1)  y'  -  t*Ay  =  taf(t),  1  <  t  <  »,  a  e  R,  a  >  -1  , 

(2.2)  ye  C( (1,-] )  , 

(2.3)  By(1)  =  8 
shall  be  approximated  by  the  'finite'  problem 

(2.4)  x,J,  -  taAxT  -  taf(t>,  1  <  t  <  T  , 

(2.5)  Bxt<1)  *  6  , 

(2.6)  SfTlx^tT)  =  y(T) 

as  T  approaches  infinity.  A  is  assumed  to  be  a  real  n  x  n  matrix  with  the  Jordan 
form  J : 

(2.7)  A  -  EJE-1  *  8 
and  J  has  the  block  diagonal  form 

(2.8)  J  -  diag(J*,j®, J-) 

where  J+  contains  the  eigenvalues  of  A  with  positive  real  part,  J°  the  eigenvalues 
of  A  with  a  zero  real  part  and  J~  the  eigenvalues  of  A  with  a  negative  real  part. 

The  dimension  of  these  three  matrices  are  r+,  r0  resp.  r_  and  the  geometrical 

multiplicity  of  the  eigenvalue  zero  will  be  called  rQ . 

The  projection  like  matrices  G+,  Gq,  G_  and  GQ  are  obtained  by  taking  the 
matrices  D+,  DQ,  D_  and  DQ,  which  are  the  projections  onto  the  direct  sums  of  invariant 

subspaces  of  J  belonging  to  eigenvalues  with  positive,  zero,  negative  real  part  resp. 

onto  the  direct  sums  of  eigenspaces  belonging  to  zero  eigenvalues  of  ,  and  by 
cancelling  all  colmuns  of  these  matrices  which  have  only  zero  entries.  So  G+  is 
n  *  r+,  Gq  is  n  x  rQ,  G_  is  n  x  r  and  G^  is  n  x  . 

By  substituting 

(2.9)  u  =  E’’y 
we  get  the  problem 

(2.10)  u'  -  t°Ju  «  taE*1f(t),  1  <  t  <  “  , 


(2.11) 


u  e  C( l 1 ,»]  > 


(2.12)  BEu(1)  «  S  . 

The  general  solution  of  (2.10),  (2.11)  is 

-  rn+r 

(2.13)  u(t)  =  [<Mt)G0,<j>(t)G_]n  +  (Hf)(t),  n  f  c  n 

where 

(2.14)  $(t)  =  exp(  — ^~1  t“+1) 
and  H  is  a  solution  operator  for  the  inhomogeneous  problem: 

Hf  =  H+f  +  HQf  +  H_f 

where 


(2.  IS) 


(2.16) 


(2.17) 


(H+f ) ( t )  =  *(t) 

( Hq  f ) ( t )  =  <p(t) 

(H_f)(t)  =  (Mt) 


/  D+$  ^sJE  'f(s)sads  , 

as 

/  Dq<P  \ s)E  1f(s)s°lds  , 

CO 

/  D_4>  ^  (  s )  E  1f(s!sads 


6 

holds  for  6  >  1.  f  is  assumed  to  be  in  C([1,“])  and  in  order  to  make  the  integral  in 
(2.16)  exist  we  assume  that 

(2.18)  «DQE'1f(t)ll  =  0(t'(a+,)r-£),  c  >  0  . 

Here  r  is  the  maximal  dimension  of  the  invariant  subspaces  of  J  associated  with 
imaginary  eigenvalues.  We  assume  that  r  >  0  because  the  case  r  *  0  has  been  treated  by 
de  Hoog  and  Weiss  (1980b).  An  analysis  of  the  operator  H  can  be  found  in  de  Hong  and 
Weiss  (1980a)  and  Lentini  and  Keller  (1980a).  Markowich  (1980b)  has  shown  the  following 
estimates,  which  hold  for  t  >  5  i  1 

(2.19)  I]  ( H  f)(t)l  <  const.  JD  E-1fl,  ,  , 


(2.20) 


M Hgf  )  ( t )  B  <  const,  t  C  maxis  1  '  r+r’ Up E  1  f  (  s  )  II  , 

s>t 


i 


(2.21) 


I (H  f)(t)l  <  const,  t  ^  max  Is^D  E  ^f(s)»  tor  y  >  0  . 
6<s<t 


All  constants  are  independent  of  f  _and  6.  Moreover  we  assume  that  B  is  a 
_  rQ+r  _  _ 

(r +  r  )  *  n  -  matrix,  that  8  e  R  and  that  the  ( rQ  +  r  )  x  ( rQ  +  r  )  matrix 


(2.22) 


BE 


($ ( 1 )G  ,♦( 1 )G_]  is  regular 


so  that  (2.3)  defines  r^  +  r  independent  boundary  conditions.  According  to  Markowich 

( 1980a, b)  these  propositions  are  necessary  and  sufficient  for  the  unique  solvability  of  the 

problems  (2.1),  (2.2)  and  (2.3)  for  all  appropriate  _f's.  Therefore  S(T)  has  to  be  a 
—  V(r0-ro> 

( r+  +  (rQ  -  rQ ) )  x  n  -  matrix  and  y(T)  f  R  so  that  (2.5)  and  (2.6)  set  up 

n  boundary  conditions. 

At  first  we  prove  a  stability  estimate  for  (2.4),  (2.5),  (2.6): 

Theorem  2.1:  We  assume  that  (2.22)  holds  and  that  (A),  (B),  (C)  which  are  defined  as 
follows,  are  fulfilled. 

(A)  1S(T)I  <  const.  for  T  ♦  "  , 

(B) 

(C) 


IS(T)EGqI  -  o(t'<01+’ )<r_1  *)  for  T 
■  -1 . 


I (S(T)EG+,S(T)EG0)  I  <  const,  for  T  ♦  » 

where  Gq  is  the  n  x  (rQ  -  rQ)  -  matrix  which  is  obtained  by  cancelling  the  columns  of 

the  matrix  D„  ■  D„  -  D„  which  have  only  zero  entries. 

0  0  0 

Then  the  problem  (2.4),  (2.5),  (2.6)  has  a  unique  solution  xT  for  all  T 

sufficiently  large  and  xT  fulfills  the  stability  estimate: 

(2.23)  •  x  ■  <  const. OB«  +  T* “+ 1 1 ( r~ 1 1 ly (T ) I  +  T<a+1)rlfl  > 

T  [ 1 ,T)  [ 1 , T ] 

Vr-  n-(Vr-’ 

if  f  e  C( ( 1 , T) ) ,  8  *  R  ,  Y(T)  e  R 

In  order  to  prove  this  we  first  reorder  J1^  by  permuting  its  lines  and  columns  so  that 


RJ°  R_1 


(2.24) 


«  3° 


-5- 


* _ aJ  « V 


where  R  is  an  appropriate  permutation  matrix*  This  corresponds  to  a  reordering  of  the 
columns  of  E*  The  reordered  matrix  will  be  called  E*  Now 


’  8  ' 

}V(Vro> 

'  6 

I- 

rn 

I  - 

ro-ro 

(2.25) 

(a)  G0  = 

0 

,_8 _ j 

t 

}r_ 

(b)  Gq  = 

U  U 

holds.  From  (2.24)  it  is  easily  concluded  that 


/4.a+1 

e/t  - 

1) 

6 

(2.26) 

exp(_4_  (t-i .  „) . 

,  a+1 
e2(t  - 

1) 

I- 

^oj 

ro  '  ro 

7o 

We  substitute 

(2.27)  vT  -  r\ 
and  get  the  problem 

(2.28)  V.J,  -  t°JvT  =  t^lT'f  (t)  ,  1  <  t  <  T  , 

(2.29)  BEv  <1 )  =  0  , 

(2.30)  S(T)Evt<T)  ■  TT (T ) 
where  J  has  the  block  structure 

(2.31)  J  =  diag(J+,J°,J~)  . 

We  write  the  general  solution  of  (2.28)  as  follows: 

(2.32)  vT(t)  =  A(t,T)C1  +  C(t)?2  +  Vp(t,T) 
where 


(2.33) 


A(t,T) 


f  , * a+1  _a+1 4  ^ 

expU  +  i  (t  ‘  T  >3 

8 

8 

-7t“+1  -  ') 

8 

e2(ta+1  -  1) 

9 

_ _ _  / 

6 

V. _  J 

J 

)  " 


and 


r„  -  r. 
0  0 


(2.34) 


C(t) 


'  8 

8 

■N 

8 

6 

, 

I 

8 

fxpla  +  1  \ 

J 

V<Vro’  r  ro+r- 

80  that  51  e  c  and  i  C 


holds,  v  (t,T)  is  an  appropriate  particular 
P 


solution  which  will  be  defined  later. 

From  (2.26)  we  easily  derive  the  following  properties  of  e^,  e2: 


(2.35) 


(a)  e(0)-I 

0  0 


(b)  e,(0)  -  e  , 


(2.36) 


(a) 


e^t01*1  -  1)  1  »  e^ ( 1  -  t“+1)  , 


(2.36) 


(b) 


e2(ta+1  -  1 ) e  1  ( 1  -  ta+1)  =  -e2(1  -  t“+1  ) 


for  all  t  e  R.  A  more  general  statement  than  (2.36)  is 
(2.37)  (a) 


ei(to+1  '  1)ei(1  -  O  “  ei(to+1  -  ti+1)  * 


a+1  a+1  a+1  a+1  a+1 

(2.37)  (b)  e2(t“  -  lie,!!  -  t“  )  =  e2<t"  -  )  -  e^l  -  ty  ) 


for  all  tg.tj  e  R. 


(2.38) 


(a) 


le(ta+1  -  1).  -  0<t<a+1,(r-,,>  for  t 


-7- 


(2.38) 


(b) 


»e2(ta+1  -  1)1  =  0(t(at1l(r-’>)  for  t  *  - 


The  estimate  (2.38)  is  derived  by  using  that 

a+1 


(2.39) 


/  k  a+1 1  1  a+1  -i  , 

e*p(^T~7  t  )  =  exp(i  t  )F(t 


where  Jk  is  an  rk  dimensional  Jordan  block  with  the  imaginary  eigenvalue  iy  and 
F(t)  is  a  real  matrix  whose  entries  are  polynomials  of  maximal  degree  ( r^  -  1)(a  +  1). 

By  inserting  (2.32)  into  the  boundary  conditions  (2.29),  (2.30)  we  get  the  linear 
block  system 


(2.40) 


fA,(T)  A  1 

c. 

8  -  bIv  ( 1 ,T) 

1 

' 

P 

A,(T)  A  (T ) 

e. 

Y  (T )  -  S (tTev  (T,T) 

.3  4  J 

P 

where 

(2.41)  (a) 
(2.41)  (b) 
(2.41)  (c) 
(2.41)  (d) 


a+1 


A1  (T)  =  [BEG^exp(  —  —  (1  -  T 

A2  =  [BEG0<BEG_exp(-  J+  — 


) ) ,BEG. 


)]  , 


~  — =  a+ 1  ~ —  a+1 

A3(T)  =  [S  (T  )  EG+ , S  (T  )EGQe  ^  (T  -  1)  +  StTlEG^lT  -  1)]  , 


~~  ~  J  a+ 1 

A4<T)  =  [S(T)EGQ,S(T)EG_exp(^— j-  T  )]  . 


-1 


The  system  (2.40)  is  soluble  iff  the  matrices  Aj  and  (Aj(T)  -  A^(T)A^  A^(T))  are 

invertible.  A^  is  invertible  by  assumption  (2.22)  and  the  existence  of 

(A,  -  A  A  1A  )  '  has  to  be  proven.  We  will  show  that 
3  4  2  1 


(2.42)  (A  -A  A-Vf1 
3  4  2  1 


0  e1  (  1  -  Ta+1  ) 


[S(T)EG+,S(T)EG0J_1(I+o(1))  for  T 


The  existence  of  the  right  hand  side  of  (2.42)  is  assured  by  proposition  (C)  of  Theorem 

2.1. 

We  split  A^(T)  into: 

(2.43) 


A3<T)  =  [S(T)EG+,S(T)EGoei(Ta+1 


1)]  +  [0,S(T )EG  e  (T 


a+1 


0  2 


-  1  >  1 


A3(T) 


Aj(T) 


and  get 


-8- 


A  (T )  =  (S(T)EG ,S(T)EG  ] 

L  9  e^T0  -  1) 


From  (C)  and  from  (2.361(a)  we  conclude 


( A3  (T  )  ) 


-!  f  1 

.  9 


( 1  -  Ta+1  ) 


|  [S(T)EG+, 


S(T)EGq) 


Aj(T)(a’(T))_1  =  [8(S(T}EGoe2(T0l+1  -  1  )e1  ( 1  -  t“+1  )  J  [S  (T  )EG+,  S  (T  )EGQ]  “ 1 

(2.46)  —  .  „  _ 

=  [9,  -  S(T)EG0e2(1  -  T  )] [S (T)EG+ ,S(T)EG  ] 

holds  because  of  (2.361(b).  The  proposition  (B)  and  (2.381(b)  assure  that 

(2.47)  lAjfTKA^fT))'1#  -  o(1)  for  T  +  »  . 

Therefore  Aj(T)”1  exists  and  can  be  written  as 

QD 

(A  (T))'1  -  (a’(T))_1(I  +  l  (-DNa^TKa'it))'1)1)  = 

i=1  3  3 

(2.48) 

=  (a’(T))_1(X  +  o(1)(  for  T  +  -  . 


(2.49)  (a) 


(2.49)  (b) 


lA  (T)l  =  o(T_<a_1 1 (r_1 ’)  for  T 

4 

•  A^(T)I  «  0(1)  for  T  +  <■> 


and  therefore 


«A/1(T)a'1A,(T)(A,(T))_1I  =  o(1)  for  T  +  » 
4  2  1  3 


because  of  (2.48),  (2.49) ,  (2*45) •  So 


(A.  -  A.a’  A,f  =  ( A  (T ) )~"  (I  +  l  (A.(T)A*  A  (T ) ( A 
3  421  3  .  .  4  2  1  3 

i=  1 


(T) )  )  ) 


=  ( A3<T)  )-1  (I  +  o(1))  for  T  +  00 

and  (2,42)  follows  immediately.  The  linear  equation  (2.40)  can  now  be  solved: 
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(2.52) 


V 

-(A  -A  A"1*  )_1A  a"1 

3  4  2  1  4  2 

5  -  BEVp(1,T) 

*2 

at 

A^Aj’VVVj’A,)"^1 

-A2lA1(A3-A4A2lAir1. 

• 

y (T)-S(T)Ev  (T,T) 

P 

Inserting  51  into  (2.33)  the  term  A(t,T)(A3(T)  -  A^(T)A~^A^(T))”^  appears.  From  (2.33) 
and  (2.42)  we  get 


(2.53)  A(t,T)(A,  -  A  a”1  A  j""1 
J  4  2  1 


(  J  ..Ct+1  Ctfl.'N 
exp(—  (t  -T  >] 

6 

6 


e  ,  (ta+1-1 )e^( 1-Ta+1 ) 
e2(ta+"-1)e1(1-Ta+1) 


0(1)  . 


Using  (2.37)  we  conclude  that 

(2.54)  I A( • #T ) ( A  (T)  -  A.(T)a‘'a,(T))''|.,  -  0 (T (“+1 H 1-1 ) ) 

3  4  21  [  i  f  T  J 

holds. 

From  (2.52),  (2.33)  and  (A)  we  derive: 

(2.55)  <  const.  (161  +  T(“+1  1  ( ^ 1  ’ ly  (T )  I  +  Iv  (1,T)I  + 

i  l'»‘J  p 

+  T(“+1)(r'1)IVp(T,T)l  +  ,VpC-T)l[1/r))  . 


Now  V  ( •  ,T)  has  to  be  defined.  We  set 
P 


(2.56) 


f  ( t,T) 


f(t),  1  <  t  <  T 

f(T),  T  <  t  <  “ 


and 


(2.57)  (TH+fXt)  t  =  /  exp(— ^  (ta+1  -  sa+1))D+E  1saf(s)ds  , 

T 

<TH0f,<t)  S”  ^  exP(~+~'  <t0+1  ~  S“+1))D0i'1sC,f(s)ds  , 


(2.58) 


(2.59)  (  H_f)(t)  -  (H_EE  f)(t),  6  <  t  <  T,  1  <  6  <  T 


so  that  we  can  define: 
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(2.60)  v  («,T)  -  Hf  -  -H  f  *  H.f  +  mH  f 

p  T  T  +  TO  T  - 

From  the  estimates  (2.19)  and  (2.21)  we  conclude  that 


(2.61)  (a) 


because 


(2.61)  (b) 


1_H  f  +  _H  fl..  <  const .  I  f  I  .  , 

T  +  T  -  [6,T]  [S,T) 


l_H  f  I . ,  <  const. IH  I , .  lf('  r)lr.  . 

T  +  [4,T]  +  [4,"1  [5,*] 


holds.  Moreover  (2.58)  can  be  estimated  as  follows 


(2.62) 


KTH0f)(t)l  <  const. lfl(t(Tj  /  (t1 


a-tl  a+1 .  ( r-1 )  a.  . 
-  s  )  s  ds  < 


<  const.  T^a+^rlfl 


[t,T]  * 


Altogether  we  get 


■H  <  c4T^°+1 >rlf I r,  , 

T  [6,T]  [o  ,T) 


(2.64)  Iv  <T,T)I  -  I (  H  f ) (T ) I  <  const. If!.,  . 

P  i  "  L'.TJ 

holds  because  (THQf)(T)  *  0.  From  these  estimates  anil  from  (2.55)  the  stability  estimate 


(2.23)  follows. 


In  order  to  derive  a  convergence  statement  we  write  the  problem  (2.1),  (2.2),  (2.3)  as 


follows 


y'  -  taAy  -  taf(t),  1  <  t  <  T,  y  «  C((1,-)) 
ByCI )  -  B  » 


(2.67)  S(T)y(T>  -  S(T)y(T) 

and  subtract  (2.4),  (2.5),  (2.6)  from  (2.65),  (2.66),  (2.67).  We  get 

(2.68)  (y  -  Kj)'  ~  taA(y  -  XT>  -  0  , 

(2.69)  B(y  -  tjHI)  -  0  , 

(2.70)  S(T)  (y  -  sty)  (T)  «  S(T)y(T)  -  7(T)  . 

If  S(T)  fulfills  the  assumptions  (A),  (B),  (C)  Theorem  2.1  can  be  applied  giving: 

(2.71)  ly  -  xTl(1  T]  <  const.  T( a+ 1  * ( r" 1 1  IS (T) y (T)  -  7(T)I  . 
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Setting  Y(T)  =  0  and  using  (2.9),  (2.13)  we  get 


(2.72)  ly  -  x  I  ,  <  const.  T( “+1 ’ ( r' 1 ’ IS (T )E (T )G„ (T )G  In  +  S (T )E ( Hf ) (T  )  I 

T  l 1 ,TJ  0  “ 


r-+ro 

for  some  n  e  c  .  Assumption  (B)  guarantees  that  all  columns  of  the  fundamental  matrix 
which  are  constant  are  dampened  by  an  o(T  *a+1^r  1 ' ) .  ah  other  appearing  columns  decay 
exponentially.  Therefore  the  term  which  originates  from  the  solution  of  the  homogeneous 
problem  converges  as  an  o(1)  as  T  *  So  Assumption  (B)  is  necessary  for  convergence 
for  general  g  and  f.  Now  let 


(2.73) 

From  the  estimates 

(2.74) 


i f ( t) i  -  o(t"(a+1)(2r_1 >_£), 


e  >  0  . 


(2.19),  (2.20),  (2.21)  we  conclude  that 

I (Hf ) (T) I  =  0(T_<a+1 )<r_1 )-E> 


holds.  Altogether  we  get 

(2.75)  ly  -  XTI(1  Tj  <  const.  (T(“+1  )(r_1  )(IS(T)EG()I  +  IIS  (T )  E<J>  (T  )G_I )  +  T_£  >  . 

If  f(t)  contains  an  exponentially  decreasing  factor  so  that  it  has  the  asymptotic  behaviour 

(2.76)  lf(t)l  =  o(t®exp(-  ^  ^  ^  t0^1)),  co  >  0 

then  there  is  an  operator  H  so  that  yp  =  EHf  is  a  particular  solution  of  (2.65)  and 

(2.77)  I  ( Hf )  (T ) I  -  o(TB+<a+1)rUnT)exp(-  q  “  <  t“+1)) 

holds.  A  proof  for  this  can  be  found  in  Markowich  (1980b).  In  this  case  T  £  in  (2.75) 
has  to  be  substituted  by  TS+<a+1  )j  exp(-  — Ta+1). 

An  optimal  choice  S(T)  =  sD(T)  would  be  so  that 

(2.78)  [Sd(T)e5q,Sd(T)EG_]  =  8 
holds.  (2.78)  is  fulfilled  for 


(2.79) 


S  (T )E  S  S  E 
D  D 


(G+ ) 


=  T 

<V 


<==>  etst  „  [G+,CJ0J 


These  are  linear  equations  for  the  rows  of  Sp  which  can  be  chosen  independently  of  T. 
The  asymptotic  boundary  condition  (2.79),  called  projection  condition  fulfills  (A),  (B), 
(C)  in  Theorem  2.1  and  is  optimal  in  the  sense  that  it  makes  the  first  two  terms  on  the 
right  hand  side  of  (2.75)  vanish. 


3.  Linear  Variable  Coefficient  Problems  -  Distinct  Eigenvalues 
In  this  chapter  we  analyze  the  problem 

(3.1)  y'  -  taA(t)y  =  taf(t),  af  nq,  1  <  t  <  »  , 

(3.2)  y  e  C<(1 ,»]  )  . 

(3.3)  By ( 1 )  =  6 
and  we  require  the  n  x  n  matrix  A(t)  to  fulfill 

(3.4)  A  «  C( [1,«] ),  A(-)  *  0  , 


(3.5)  Alt)  *  \  A.t  1  for  t  sufficiently  large  . 

i»0  1 

Moreover  let  JQ  be  the  Jordan  form  of  Ag  obtained  by 

(3.6)  Aq  -  EJ0E-1  . 

The  following  assumption  is  basic  for  this  chapter: 

(3.7)  =  diag( X, » • . .  ,X  ),  X.  *  X.  for  i  *  j  . 

o  1  n  i  j 

The  substitution 

(3.8)  u  =  E_1y 
gives  the  problem 

(3.9)  u'  -  t°\j(t)u  -  taE_1f(t),  1  t  t  <  «■  , 

(3.10)  u  e  C( [1,»J ) 
where 

OB 

(3.11)  J(t)  =  E~ 1  A(  t )  E  =»  l  J,  -  E~Ve  . 

i“0  1  1  1 

The  fundamental  matrix  of  the  homogeneous  problem  (3.9)  can  be  represented  as  an  asymptotic 
series  (see  Wasow  (1965)  and  Coddington  and  Levinson  (1955)): 

(3.12)  $(t)  =  P(t)t°e2(t) 
where 
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(3.13) 


P(t)  ~  I  +  l  P  t*1  , 
i-1 

(^•14)  D  =  diag( dj, . . . ,dn)  , 

t“+1  ta 

(3*'5>  Q(t)  “  JQ  ”-+  -  +  Q1  —  +  •••  +  Qat»  are  diagonal  matrices 

hold.  The  unknown  coefficients  P^  and  D  can  be  calculated  by  algebraic  operation 
from  the  J^s.  An  algorithm  for  that  is  given  in  Markowich  (1980b),  and  therefore  the 
asymptotic  behaviour  of  the  basic  solution  can  be  determined  knowing  the  Q.' s  and  D. 

Let  D0  ^  the  Pro3ection  onto  the  direct  sum  of  eigenspaces  of  JQ  associated  with 
those  eigenvalues  with  a  real  part  zero  which  produce  a  basic  solution  which  is  in 
C< [1,-1 )  and  let  GQ  be  the  projection  like  matrix  which  is  obtained  by  cancelling  those 
columns  of  DQ  which  have  only  zero  entries.  Gfl  be  a  n  x  rQ  matrix.  Then  the  general 
solution  of  the  problem  (3.9),  (3.10)  is 

~ 

(3.16)  u(t)  =  [<Mt)G0,<Mt)G_]0  +  (Hf)(t),  n  c  c 

where  G_  is  defined  as  in  Chapter  2  and  up(t)  “  (Hf)(t)  is  an  appropriate  particular 
solution,  which  has  been  described  by  Markowich  (1980b).  The  operator  H  operates  on  the 
space  of  all  functions  fulfilling 

(3.17)  f  e  C([6,<»]),  S  >  1  and  lf(t)l  =  Ott'01"1"6),  e  >  0  . 

Then  the  estimate 

(3.18)  1 ( Hf ) ( t ) I  <  const,  t  Efnt  •  maxisa+1 +E f ( s ) I 

s>6 

has  been  proven  by  Markowich  (1980b).  The  particular  solution  on  (I,”)  is  obtained  by 

continuation.  The  boundary  value  problems  (3.1),  (3.2),  (3.3)  is  -  under  the  given 

r  +r 

assumption  on  Aft)  and  f  -  for  all  B  6  R  0  “  uniquely  soluble  iff  the 

(rQ  +  r_)  x  (rQ  +  r_)  matrix 

(3.19)  BE  [$  ( 1  )Gg  ,<p  ( 1  )G  ]  is  regular  . 
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Of  course,  B  has  to  be  a  ( rQ  +  r  )  x  n  matrix.  We  consider  the  approximating  problems 

(3.20)  x^,  -  t*A(t)xT  =  t“f(t),  1  <  t  <  T,  a  e  nq  , 

(3.21)  Bxt(1 )  -  B  , 


(3.22)  S(T)xt(T)  -  y(T) 

^  1 
S (T )  is  a  (n  -  (rQ  +  r_)  >  x  n  matrix  and  -y(T)  €  R 


n-(r0+r_) 


For  the  following  G+  is 


defined  as  in  Chapter  2  and  GQ  is  the  n  x  (rQ  -  rQ)  matrix  which  is  obtained  by 
cancelling  the  zero  columns  of  -  D^.  Then  the  following  stability  theorem,  which  is 
the  analogous  to  Theorem  2.1,  holds : 

Theorem  3.1;  Assume  that  (3.19)  and  <A1(  ( B ^ ) ,  (C^)  which  are  defined  as  follows,  holds 
(A.)  IS(T)I  <  con6t.  as  T  +  •  , 


IS(T)EGq I  =o(1!  as  T  ♦  “  , 

m  —  1 

I [S(T)EG+,S(T)EGq]  I  5  const,  as  T  ♦  “  . 


Then  the  problems  (3.19),  (3.21),  (3.22)  has  a  unique  solution  xT  for  sufficiently 
large  T.  x_  fulfills  the  estimate 


(3.23) 


lx  I  <  const. (161  +  ly(T)l  +  T  inTlfl.  .) 

T  1>»TJ  llf*J 


vr-  n-<r0+r-> 

for  f  e  C([1,T]),  B  e  R  ,  1 (T )  e  R 


The  substitution 


gives  the  new  problem 


v^,  -  t“j(t)vT  =  taE_1f(t),  1  <  t  <  T  , 
BEvT(1)  «  B  , 


(3.27)  S(T)EVT(T)  -  YtT)  . 

As  the  general  solution  of  (3.25)  we  take  for  convenience 

(3.28)  vT(t)  =  $(t)e"2!T)T'DlG+,G0JC1  +  1(t)  (G0,Gjt2  +  v  (t,T) 


r++(r0-r0)  r0+r- 

where  e  C  ,  e  C  hold  and  v  (*,T)  is  an  appropriate  particular 

1  2  p 

solution  which  will  be  defined  later.  <)>  ( t )  is  the  fundamental  matrix  as  of  (3.12). 
Evaluation  of  the  boundary  conditions  (3.26),  (3.27)  gives  the  linear  block  system 


(3.29) 


‘BE4>(1)e"C(T)T-D[G+,G0] 

BE*(1) [G0,G_] 

[V 

’  8  -  BEv  (1,T) 

P 

_  S(T)EP(T)(G+,G0] 

S(T)E*(T) [G0,G_J_ 

s2J 

__ 

Y (T)  -  S(T)Ev  (T,T) 

P 

The  matrix  in  the  (1.1)  position  is  bounded  because  of  the  definition  of  G+,  Gq  and 
because  of  the  diagonal  form  of  Q(T)  and  D.  The  matrix  in  the  (2.2)  position  is  bounded 
too,  because  $(t)[GQ,G  ]  is  the  matrix,  whose  columns  are  the  basic  solution  of  the 
homogeneous  problem  which  are  in  C([1,“J)  and  because  ( A ^ )  holds.  Moreover 

(3.30)  IS(T)E<(>  (T)  I  =  0(1)  as  T  ♦  “ 

because  of  (B^)  and  (3.12).  The  matrix  in  the  (1.2)  position  is  invertible  because  of 
(3.19)  and  its  inverse  is,  as  the  matrix,  independent  of  T.  Finally 

(3.31)  S(T)EP(T) [G+,Gq 

because  of  the  asymptotic  expansion  for  P(t),  (C.j)  assures  the  bounded  invertibility  of 

the  matrix  in  the  (2.1)  position. 

From  (2.40)  and  (2.52)  we  conclude  immediately  that  the  system  (3.29)  has  a  unique 


S (T )E [G+ ,Gq ]  +  0 (T~ 1 ) 


solution  (€.|/£2)  €  C  and  the  estimate 


(3.32) 

follows. 


,Vt'(1,T]  <  con8t'(,e11  +  +  ,vpC-T)l[1(T]  ) 


The  particular  solution  v^(»,T)  has  to  be  defined  now.  We  set 


(3.33)  (a)  v  C,T)  -  H^’f  -  H^’f  ♦  H^’f  +  H(1)f,  (b)  h! '  ’  f  =  £  Ji’11! 

p  T  T  +  TO  T  -  TO  .  ,  T  0 1 

1=1 


where 


(3.34) 


(^j^fHt)  =  <J> ( t)  /  D+i)i'1(s)saE'1f(s)ds  , 


-1...  _a_-1. 


(3.35)  (  H(1)f)(t)  =  <Mt)  /  D  <T'(s)s'V1f(s)ds,  «  >  1 

T  "  5 


and 
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if  (I)  holds 


(3.36) 


<THdrf)<t) 


$(t)  /  1(s)s°l£  Vfsjds, 

T 

4>(t)  /  Dg^<)>  '(s)saE  ^(sfds,  if  (II)  holds  . 


DQi  is  the  projection  onto  that  eigenspace  of  Jg  which  belongs  to  the  ith  eigenvalue 


with  real 

part  zero  and  ( I ) 

,  (II)  are 

defined  as 

(I) 

Re(Q(t)D0i) 

+  +“  or 

(Re(Q(t)DQi)  5  0  and 

Re(DD0i)  > 

0)  , 

(II) 

Re(Q(t)D0i) 

♦  or 

(RE(Q(t)DQi)  5  0  and 

Re(DDQi)  < 

0)  . 

From  the 

considerations  in 

Markowich 

(1980b),  Chapter  3,  we 

immediately 

conclude  that 

(3.37) 

IH(1 >f« 

.  ,  ,  <  T“+1tnTlf 1 , .  , 

T 

(6  ,T]  (6  ,T) 

Therefore  the  estimate  (3.23)  follows  and  Theorem  3.1  is  proven. 


As  in  Chapter  2  the  convergence  estimate  follows 

(3.38)  «y  -  <  const. IIS(T)y(T)  -  y(T)« 

for  all  (r+  +  (rQ  -  rQ )  )  *  n  matrices  S(T)  which  fulfill  (A.,),  ( B  ^ )  and  ( C 1 ) . 

Setting  y(T)  =  0  and  inserting  (3.16)  we  conclude 

(3.39)  ly  -  XTl;i  T]  <  const.(IS(T)E*(T)[G0,G_JI  +  IS (T)E(Hf ) (T) I >  . 

The  assumptions  ( A ^ )  and  ( B 1 )  guarantee  convergence  for  all  f  fulfilling  (3.17)  because 
(3.18)  holds.  If 

(3.40)  S(T)EGq  =  0 

convergence  of  the  order  T  ^  +  T  E£nT  follows.  In  many  practical  cases  all  eigenvalues 
with  real  part  zero  produce  exponentially  decaying  solutions  and  f  also  decays 
exponentially.  The  operator  H  can  be  changed  to  an  operator  H,  so  that  (Hf)(t) 
decays  with  the  same  exponential  factor  (see  Markowich  (1980b)).  In  this  case  exponential 
convergence  follows  from  (3.39). 

The  optimal  boundary  condition  is  again  the  projection  condition  and  it  has  to  be 
calculated  from  the  equation 
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(3.41) 


Sd(T)£P(T) 


fG 


T 

+ 


‘V 


which  is  uniquely  soluble  because  of  the  regularity  oj.  E  and  P(T). 

The  asymptotic  boundary  condition 

(3.42)  SjjfTlJtj.fT)  »  0 
would  imply  that 

(3.43)  Sd<T)E$(T)  [G  ,GJ  S  e 

because  of  the  form  of  4>(T).  However,  we  do  not  know  P(T),  but  we  can  calculate  the 
coefficients  P^  of  its  expansion  recursively  (for  the  algorithm  see  Markowich  (1980b)). 
Having  calculated  P^ ,P2« • • • ,Pk  we  set 


(3.44) 


P(T) 


I  + 


k 

I 


i-1 


PiT 


i 


and  solve 

(3.45)  iQ(T)EP(T)  » 

instead  of  (3.41).  Because 

—  -k-1 

(3.46)  PIT)  »  P (T)  +  0(T  )  as  T  ♦  ■ 

holds  we  get  by  a  simple  perturbation  analysis 

(3.47)  IQ(T)  *  sd(t>  +  as  T  +  »  . 

Therefore 

(3.48)  ISD(T)E<MT)  [G0,G_]  I  <  const.  t"*"1 1$ (T )  [GQ ,G_)  I 

holds  and  the  boundary  condition  SD,T*XT*T)  =  0  implys  at  least  convergence  of  the 
order  T-Ic_1  if  f  =  0  holds.  More  generally  speaking  the  order  of  convergence  is 
determinied  by  inserting  (3.48)  into  (3.39). 

However,  this  rather  work-intensive  procedure  does  only  make  sense  if  some  columns  of 
$(t)[Gg,G  )  do  not  converge  exponentially.  Only  in  this  case  the  projection  conditions 
imply  a  significant  improvement  of  the  order  of  convergence. 
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4.  Linear  Problems  -  The  General  Case 


In  this  chapter  we  admit  a  general  Jordan  form  of  AQ.  So  we  deal  with  problems  of 
the  form  (3.1),  (3.2),  (3.3)  with  the  assumption  (3.4),  (3.5),  (3.6). 

Again  we  perform  the  substitution  (3.8)  and  get  (3.9),  (3.10). 

The  fundamental  matrix  <f>(t)  of  the  homogeneous  problem  (3.9)  can  now  be  represented 
as  an  asymptotic  log-exponential  power  series  of  the  following  form 
(4.1)  4>(  t)  =  P(t)tDeC(t) 

where 


(4.2) 


p<t)  ~  i 


i=0 


P.  t 
1 


P  £ 


(4.3) 


is  a  constant  matrix  in  Jordan  form 


(4.4) 


“+1-  ;  i 

Q(t)  -  diag(J0>  +  Ql  ^1).1  +  •••  +  Cp(a+1).,tp 


and  the  are  diagonal  matrices.  diag(J0)  is  the  matrix  which  has  the  same  diagonal 

entries  as  JQ  and  all  other  entries  zero.  The  matrices  T°  and  e^*1*  commute  because 
the  diagonal  elements  of  Q(t)  which  belong  to  a  particular  Jordan  block  of  D  are 
equal.  Moreover  P(t)  can  be  split  up  into: 


(4.5) 

where 


Pit) 


P(1)(t) 


P(2»tt> 


P( 3 )  < t) 


(4.6) 


(1)(t)  ~  1  +  J.  P(1)it'i  * 


i-  1 


(4,7) 


P(k)(t)  ~  L  P(k)it  '  *  -  2'3 

1=0 


P^2)(t),  P^j(t)  are  *n  bl°ck  diagonal  form  too.  The  i-th  diagonal  block  of  P(2>(t) 
corresponds  to  that  block  in  Jq  which  is  obtained  by  gathering  all  Jordan  blocks 
belonging  to  the  i-th  eigenvalue  of  JQ  and  the  j-th  diagonal  block  of  P ( 3 ) < t ) 
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corresponds  to  the  j-th  eigenvalue  of  Q(t)  where  in  both  cases  only  different  eigenvalues 


are  counted.  Markowich  (1980b)  has  shown  that 

(a+1 )  ( r  -1 ) 

(4.8)  ,<P(2)(t>)  Di*  <  const,  t 

holds,  where  is  the  projection  onto  the  direct  sum  of  invariant  subspaces  associated 

with  the  i-th  eigenvalue  of  Jg  and  r^  is  the  algebraic  multiplicity  of  that  i-th 
eigenvalue.  The  statement  (4.8)  holds  for  the  matrix  p(2)  derived  as  in  Markowich 
(1980b) . 

The  matrix  Pjjj(t)tD  is  the  fundamental  matrix  of  the  system 

(4.9)  z'  ”  (~  B  +  ~  B(x))z,  B  e  C(  [1,“]  ) 

X 

where 

(4.10)  u  -  P(1 )(t)P(2)(t)ec(t)z 

has  been  set.  (4.9)  has  a  singularity  of  the  first  kind  of  t  «  °°.  obviously  Pjjj(t) 
and  D  are  not  uniquely  defined,  only  their  product  is  unique  (neglecting  multiplication 
with  a  constant  matrix  from  the  right  side).  Pj3j(t)t°  =  (P^ltjt  ^t0*1  would  also  be  a 
way  of  splitting  the  product.  The  algorithm  given  by  Wasow  (1965)  establishes  a  matrix 
P(3)  (t)  which  has  a  convergent  power  series  expansion,  but  P^jt”)  is  not  regular.  We 
will  show  now  that  a  representation  can  be  given,  so  that  P  j  («>)  is  regular.  Therefore 
we  assume  that  B  is  in  Jordan  -  canonical  form: 

(4.11)  B  =  diag(B1,...,Bq) 

and  B^  has  the  only  eigenvalue  b^,  where  Refb^)  4  0  for  1  4  i  4  s  and  Re(b^)  >  0 
for  s  +  1  4  j  4  n.  We  write  (4.9)  as 

(4.12)  z'  =*  -  Bz  +  -  (B(x)z),  B(x)  =  ~  B(  x ) 

XX  X 

and  set  for  1  4  i  4  s 


(4.13) 


zt  ( t ) 


9  1 


9 


9 


+  (GBz  >(t>  , 


5  4  t  <  ” 


9 


where  G  is  the  operator  defined  in  Markowich  (1980b),  Chapter  4  which  applied  to  a 
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I 


(4.22) 

and  (4.12)  becomes 

(4.23) 

Now  we  set 


V1 

z .  »  z  •  t  J 
3  3 


z'  »  —  (B  -  (b,  +  1)I)z.  +  —  B(x)z 
J  *  3  3  x  j 


(4.24) 


z^(t) 


B  -<b  +1 )I 
t  3  3 

-3 

0 


+  (GBz  )(t> 


so  that 
(4.25) 

and  (4.13),  (4.21)  implies  that 


o  »  1 1 
3 


B ,-(b.+1 >1 
3  3 


I  =  c*  t  1 ( Jnt) 


dim( b . ) — 1 


(4.26) 


2  j+dim(  B  . ) 

z.(t)  -  (I  +  0(t  <lnt)  3  )) 


B  -(b  +1 >1 
t  3  3 

0 
• 

0 


holds  and  from  (4.22)  we  conclude 


(4.27) 


2j+dira(B.) 

z.(t)  *  (I  +  0(t  (int)  3 

3 


f  9 
• 

0 

B  . 


)) 


as  t  ♦  *> 


Obviously  the  matrix 

(4.28) 

is  a  fundamental  matrix  of  the  system  (4.9).  Therefore  p ( 3 ) < t )  and  D  in  (4.5)  and 
(4.1)  can  be  chosen  so  that 

(4.29) 
holds. 


z(t)  *  [z(t),...,z(t)]  -  (I  +  0( t”1 ( Znt )m) ) tB,  me  N 


IP  ) (t)  1 1  =  II  +  o(1)l  <  const. 
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Knowing  the  fundamental  matrix  asymptotically  we  can  sort  out  the  basic  solution  <t>. 
fulfilling  ».  f  C([1,“])  so  that  the  general  solution  of  (3.90,  (3.10)  is 

Vr- 

(4.30)  u(t)  =■  *<t){G  ,G_]n  +  (HfHt),  h  e  C 

where  H  is  defining  an  appropriate  particular  solution  Hf  on  [$,“)  if 

(4.31)  f(t>  =  t'(a+1,r"eUnt)iF(t),  F  e  C  (  [5  ,») ) 

b 

where  r  is  the  maximal  algebraic  multiplicity  of  eigenvalues  of  which  have  real  part 


zero.  Moreover 


j  +1 

(4.32)  I  ( Hf )  ( t )  B  *  const.  t~6Unt)  °  U  FB  .  »  ,  #  0  <  H  4  n  . 

[<5,»J  0 

The  particular  solution  on  [1,“]  is  obtained  by  continuation. 


A^»-.n  the  boundary  value  problem  on  the  infinite  interval  is  uniquely  soluble  for  all 
r  *r  ~  __ 

B  r  R  and  f‘s  which  fulfill  (4.31)  iff  the  (rQ  +  r  )  x  (r  +  r  )  matrix 

(4.33)  BE[$(  1  )GQ  ,<M  1  )G  )  is  regular  . 

Of  course  B  is  an  (r^  *  r  )  *  n  matrix. 

The  approximating  problems  have  the  form  (3.20),  (3.21),  (3.22).  We  will  again  prove 
a  stability  theorem. 

Theorem  4.1.  Assume  that  (4.33)  and  (A2),  (B^),  (C2)  which  are  defined  as  follows,  hold: 

(A2)  US(T)II  4  const,  as  T  ♦  «•  , 

<B2)  »S(T)EP(T)G0«  =  o(T  '  )  , 

(C2)  MS(T)EP(1  }  (T )  [G+,P(2)(T)P(3  (T)G0]  )_1  I  =  0(T(Q+1  ’  tr_1  f)  . 

Then  there  is  a  unique  solution  x,p  of  the  problem  (3.20),  (3.21)^  (3.22)  for  T 

r0+r_ 

sufficiently  large  and  the  following  estimate  holds  for  all  0  £  R  , 

(VV+I> 

T(T)  f  R  ,  f  t  C(  [  1 ,T]  ): 

(4.34)  «XT«(1  T]  <  const.  (IIB#  +  T(a+1  )!r-1  ’  lly(T)ll  +  T(<5+1  ’  r  UnT )  j»  f »  ( ;  T)  )  . 

We  substitute 


xT(t)  =  EP(1)(t)wT(t),  wT  = 


and  get  three  separate  problems 


+ 

1 

+ 

J  (t) 

w  ( t ) 

T 

(4.36) 

«T(t) 

=  t° 

Cl 

o 

rt 

w°(t) 

T 

WT(t) 

j"tt) 

w  (  t) 

T 

+  tap"!  (t)E_1f!t) 

V  1  _ _ 

f(t) 


where 


(4.37) 


J+(“) 


J°  (“) 


J  (») 


and  the  eigenvalues  of  J*  have  positive  real  part,  the  eigenvalues  of  have  a  zero 

real  part  and  the  eigenvalues  of  have  a  negative  real  part.  This  structure  can  always 

be  obtained  by  reordering  the  columns  of  E.  Now  we  rewrite  the  equation  for  w*: 


(4.38) 


w*(t)’  =  t“jgW*(t)  +  (J+(t)  -  J*)w*(t)  +  t“f+(t) 


We  define  the  general  solution  of  (4.38)  as: 

(4.39)  w*(t>  =  exp(^-2-j-  (t“+1  -  t“+1))5+  +  <tH+(J+  -  jJ^Ht)  +  ( t) 


0  + 

where  TH+  is  defined  in  (2.57)with  E  =  I  and  J  =  JQ. 


We  derive 


w*(t)  =  *+«t,T)?+  +  i|»+(f+)(t)  e  C([S,T]) 


(4.40)  ((I  *tH+(J+  -  J*))w*)(t)  =  exp(^-py  (t“+1  -  t“+1))C+  +  (TH+f+)(t)) 

de  Hoog  and  Weiss  (1980b)  have  shewn  that  (I  -  TH+(J+  -  J* ) )  is  invertible  as  operator 
C([5,TJ)  with  6  and  T  sufficiently  large,  so  that 

(4.41 ) 
where 

(4.42) 
and 

(4.43) 

Moreover,  they  have  shown  that 


*+C,T)  -  (I  -  T«+(J+  -  jQ))_1exp(^-pY  (h  -  Ta+1)),  h(t>  =  t 


TVV  “  (I  -  TVJ+  “ 


a+1 


-24- 


(4.44)  (a)  ‘V'TJ'ta,,.)  <  const.,  (b)  V+(  V  *  (6  ,Tl  ‘  C°nSt* *f  ‘  [6 ,T1 

and  from  (2.57)  and  (4.42)  we  derive  that 

(4.45)  i|<+(T(T)  =  I  . 

Now  we  define  the  general  solution  of  (4.36)  as 


«T<t) 


■ 

'  <^+(t,T) 

e 

G  ,P„1(t)P,„(t)eC(t,-«(T)a\ 

. 

9 

e  , 

+  (2)  (3)  0 

(4.46) 


P(2)U)P(3)(t)e0(t>tDcVG-,52  + 


51  + 


rvv(t) 


w°(t,T) 

P 


w  (t,T) 
P 


V  c 


WV 


52  «  c 


Vr- 


where  w^ ,  w^  are  appropriate  particular  solutions.  This  solution  is  defined  on  [S,T] 
and  the  corresponding  solution  on  [1,T]  is  obtained  by  continuing  iji+(*,T),  Resubstituting 
in  (4.35)  and  evaluating  at  the  boundaries  sets  up  the  linear  block  system  for  5  ,  : 


(4.47) 


BEP(1,(1) 


[V1- 


T)  9 
0  0 


VP(2,<1)P(3,<1)e«(1>-C(TVDSc 


S(T)EP1 (T) (G+,P(2) (T)P(3) (T)GqJ 


BE4>(1  )  [G0,G_] 


s(t)e*(t)  [gq,gj 


8  -  BEP (  ( 1 ) 


TW(1) 


w  (1,T) 
P 


w  (1,T) 
P 


Y(T)  -  S(T)EP(1)(T) 


r^+(f+)(T) 


w°(T,T) 

P 


* 


w  (T,T) 


The  matrix  in  the  (1.1)  position  is  bounded,  BE$( 1 ) [GQ,G_]  is  independent  of  T  and 
invertible  because  of  (4.33),  the  matrix  in  the  (2.1)  position  fulfills  (C2)  and 

(4.48)  I S(T )£$ (T ) (Gg ,G_ 1  I  <  o( 1 )T_<“+1 1 ( r_1 ’ 
holds  because  of  (Bj).  From  (2.40),  (2.52)  we  conclude  that 

(4.49)  lx  I  ,  <  const. (181  +  T^°+1^r  ^ I Y (T) I  +  Iv  (1,T)I  + 

T  l i #TJ  p 

+  t(01+1  1  (r_1  1  IS(T)Ev  (T,T)  I  +  Iv  (  •  ,T)  I  ,,  ) 

P  p  M,T) 


v  ( t,T)  i=  P... (t)  w  ( t,T) 
P  ( 1 )  P 


(Hj2)f)(t) 


I  w  (t,T)  J 
1  P 

has  been  used.  Splitting  H^2  *  into  TH|2^'  TH *2  * ,  TH^2*  where  TH^.2  *  is  already 
defined  by  _t|i  (f  ),  we  can  define  h!2  *  as  we  have  defined  in  (3.36)  only  the 

Dq  j  have  to  be  substituted  by  the  projections  onto  the  invariant  subspaces  of  D.  The 


(4.51)  'TH02)f,(1,Tl  <  T(0l+1)r(fnT)  °lfl[1jT] 

results  as  in  chapter  3. 

(2) 

tH_  f  can  be  constructed  by  the  same  perturbation  approach  we  used  for  the 
(2) 

construction  of  f .  We  set 

T  + 

(4.52)  T*J?J  =  (I  -  -  VJ'V.f.  =  w*C,T) 

where  is  defined  in  (2.59)  with  E  =  I  and  J  =  J  and  then 


(th;  fxt>  =  p(1)(t)  6 

TMt.)U) 


holds.  Moreover  the  estimate 


12  ) 

I  Hl  'fl ,,  <  const . I f I 

T  -  [  1  ,T]  !  1  ,T) 


is  fulfilled.  Because  of  (Bj)  and  (3.36)  we  get 


(4.55) 


IS(T)E<THg2)f > (T )  I  <  const.  t“+1 (4nT)  °»fl ^  T) 
so  that  the  estimate  (4.34)  follows. 

Again  the  convergence  estimate  follows: 

(4.56)  Jy  -  x^.1^  T]  <  const.  T<a+1 1  * r_1  'iSfDyfT)  -  y(T)I 

for  all  matrices  S(T)  fulfilling  (Aj),  (Bj).  (C^).  Setting  Y(T)  =  0  and  using  (4.30) 
we  get  the  order  of  convergence  as  follows: 

(4.57)  ly  -  >^l(1  T]  <  const.  T(a+1 > 1 r_1 > ( |S(T)E$ (T) (GQ ,GJ I  +  l(Hf)(T)l)  . 
Assumption  (C2)  guarantees  convergence  for  all  f'  which  fulfill 

(4.58)  lf(t)l  -  0(t"(a+1 > (2r“1 )_C),  E  >  0 


because  the  columns  in  $(T)  which  may  be  constant  as  T  ♦  “  are  dampened  by  the  factor 
o(T  Again  if  all  columns  of  $(t)  and  f  decay  exponentially,  the 

convergence  is  exponential,  too. 

Still  the  question  has  to  be  answered  whether  there  is  a  matrix  S(T)  fulfilling  the 
assumption  of  Theorem  4.1  and  how  it  can  be  constructed.  We  set 

(4.59)  S(T)EP(1)(T)  •  i(T) 
and  choose  SIT)  so  that: 

,T 

I 

(4.60) 


S(T )  [G+  ,9,0] 


e 


(4.61) 

and 


s(T)[e,e,G  ]  -  e 


(4.62 ) 


S(T)P(2)(T)P(3)(T)10,Go,0]  »T 


(a+ 


i )  ( r— i )  r  e 


’  e  ■ 

L  <«n»  J 


because  of  the  block  structure  of  p(2)>  p(3)  (4.62)  is  equivalent  to 

0 


(4.63 ) 


S(T )  (0  ,GQ,  8]  =  T 


•  ( a+ 1 )  ( r- 1 ) 


•  T  . 

<V  -> 


PU)(T)P(2)(T)  * 
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The  equations  (4.60),  (4.61),  (4.63)  determine  s(T)  and  s(T)  can  be  calculated  from 

(4.64)  S(T)  -  SfTJp’JjdJE-1  . 

SIT)  fulfills  (B2)  because  of  (4.62),  the  proposition  (A2)  follows  from  (4.8)  and  (4.29), 
and  (C2)  is  implied  by  (4.60),  (4.62).  This  asymptotic  boundary  condition  (with 
Y(T)  =0)  is  the  projection  condition  fulfilling 

(4.65)  SD(T)E(MTHG0  ,GJ  -  0  . 

In  general  we  only  can  determine  a  finite  number  of  coefficients  of  the  expansion  of  p(i)» 
P(2)»  p(3)-  An  algorithm  is  given  in  Markowich  (1980b)  and  it  is  shown  that 


(4.66) 


P(2,(t) 


M  s.(t)EipJ2)(t) 


where  the  matrices  are  in  block  diagonal  form  and  their  diagonal  blocks  are 


(4.67)  Si^(t)  «  diag(1,t  9i0,...,t  *3  9i3)  . 

The  E^  s  are  regular  and 


(4.68) 


I  + 


00 


I 

3-1 


as 


t  00 


holds • 


We  denote  by 


(4.69) 


P(2,0(t) 


m 

i  i 

i-i 


Si0(t)EiP(2)0(t)' 


Pd)0(t> 


3-1 


P13  t 
(2)0 


i 

p 


that  diagonal  block  of  p(2)  which  is  associated  with  the  zero-real  part  eigenvalue  of 

i,  ik2 

Jq.  Assume  now  that  we  know  P(2 )o '  " 1 ,P(2 )C  ^or  ^  and 

k2  >  2p(a  +  1)(r  -  1)  -  1  then  we  set 


(4.70) 


P(2)0(t) 


i-1 


**iO^^EiP(2  )0(  ^ ' 


F1 

(2)0 


P1J  t 
(2)0 


.  2 
P 


il 
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and  assume  that  we  know 


(4.71) 


P(3)0(t> 


I 


k3  >  p(a  +  1 ) ( r  -  1) 


1  . 


P(3 ID (t)  consists  of  the  first  k3  +  1  summands  of  P ( 3 )o ( D  which  is  th«t  block  of 
p^3j(t)  associated  with  the  real  par^  zero  eigenvalues  of  JQ.  Moreover  we  have  to  know 


(4.72) 


Pd,(t) 


I  + 


1 

V 

U 

j=l 


P(1)jt 


k  >  (a  +  1 ) ( r  -  1) 


1  . 


Using  p(1)'P(2)o'P(3)0  instead  of  p( 1 ) < p( 2 )0 ' p( 3 )0  we  emulate  ®  matrix  SIT)  instead 
of  S(T)  using  (4.60),  (4.61),  (4.63),  (4.64)  and  a  perturbation  analysis  shows  that 
(4.73)  I(T)  =  SIT)  +  0(T'k_1) 

where 


k  +  1  _  k  +  1 

(4.74)  k  -  min(  k  ,  — - (<X+  1)(r  -  1)  -  1,  - - - 1)  >  <ct  +  1 )  ( r  -  1) 

Ip  P 

Therefore 

(4.75) 


-k-1 , 


S(T)((i(T)  [GQ,G_]  =  0(T  )*(T)[G0,G_) 


(a+1 ) ( r— 1 ) — k— 1 


and  the  order  of  convergence  for  homogeneous  problems  is  at  least  T 

The  requirement  that  A(t)  is  analytical  in  t  =  “  is  very  restricting,  therefore  we 
will  now  admit  matrices  A(t)  fulfilling 

(4.76)  A<  C(  [1,-J  ),  a(~)  £  Clat111+,([0,  j]),  6>1 

where  l  is  the  maximal  algebraic  multiplicity  of  an  eigenvalue  of  A(")  with  nonpositive 


real  part.  Therefore  A  can  be  expanded: 


(4.77) 


A(t)  -  AQ  +  t"1A1  + 


+  t’(a+1,lVi)t  +  *(t) 


where 

(4.78)  A(t)  -  a<  t)t-<a+1  e  >  0,  af  Cfe  ( [  1  ,«■>)  >  . 

The  problem  (3.1),  (3.2)  can  now  be  rewritten  as 
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(4.79) 


y*  “  tB  f  At  'y  -  taA(t)y  +  taf(t)  , 
i-0 

(4.80)  y  e  C(  (1,“1) 
and  can  be  regarded  as  a  perturbed  system  of 

(a+1)7 

(4.81)  y'  -  t  l  A . t  y  «  taf(t)  , 

i-0  1 

(4.82)  y  «  C( II )  . 

«v  Q 

Markovich  (1980b)  has  proven  that  the  n  x  ( r^  +  r_)  solution  matrix  Et|/  of  (4.79), 
(4.80)  fulfills  for  large  t: 

(4.83)  l*°(t)  -  $(t) [Gq,G_) I  <  const,  t-1~8 ( £nt ) j 1$ ( t) [G0 ,G_] I 

where  0  <  j  <  2n  holds  and  4>(tHGQ,G_l  is  the  general  solution  of  (4.81),  (4.82). 
Moreover  a  particular  solution  E(i(f)  of  (4.79),  (4.80)  can  be  constructed  if  f  fulfills 
(4.31)  and 

(4.84)  li(r(f)(t)  -  (HfXt)l  C  const,  t”1"8  (tnt) j  1  (Hf )  ( t )  I 

r.+r_ 

holds.  The  problem  (3.1),  (3.2),  (3.3)  is  for  all  6  e  R  and  f  fulfilling  (4.31) 

uniquely  soluble  iff  the  ( rQ  +  r_)  »  (rQ  +  r_)  matrix 

(4.85)  BE<(i_(1)  is  regular  . 

Of  course,  B  is  a  (r^  +  rjxn  matrix. 

For  the  existence  theory  and  for  the  following  stability  theorem  it  is  sufficient  to 
require  that  (4.77)  holds  with  £  *  r  and  with  Aft)  *  t  Ca(t).  This  implies  the 

right  hand  sides  of  (4.83)  and  (4.84)  to  equal  const .( inT ) ^T  6. 

We  consider  the  asymptotic  boundary  value  problem 

(4.86)  x^,  -  t°A(t)xT  «  taf(t),  1  <  t  <  T  , 

(4.87)  Bx  (1 )  =  S  , 

(4.88)  S(T)>^(T)  »  Y(T) 

and  show  that  the  construction  of  S(T)  and  the  stability  estimte  (4.34)  depend  only  on 
the  validity  of  (Aj),  ( > »  (Cj)  for  the  perturbed  problem 
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(4.90) 


S  (T )  x^,  (T )  -  Y(T) 


if  (4.85)  holds. 

Theorem  4.2.  Let  A  fulfill  (4.76)  and  let  (4.85),  (Aj),  (Bj),  (Cj)  of  Theorem  4.1  hold 
where 

*(t)  -  P(1)(t)P(2)(t)P(3)(t)tReS(t) 

is  the  fundamental  matrix  of  the  homogeneous  problem  (4.81).  Then  there  is  a  unique 

solution  xT  of  the  problem^ ( 4 . 86 ) ,  (4.87),  (4.88)  for  T  sufficiently  large  and  for  all 
r  +r  r  +(r  -r  ) 

6  e  R  ,  Y  (T )  e  R  ,  f  «  C((1,T]).  This  solution  x,j,  fulfills  the  estimate 

(4.34). 

He  write  (4.87) 


(4.91  ) 


(a+1 ) r 

*(  l 

i-0 


Ait 


Cl  ~ 

t  (A(t)x 


+  f(t)). 


1  <  t  <  T 


and  write  the  general  solution  after  having  set  xT  -  EvT  as 


(4.92) 


vT(t)  -  P, 


(i) (t>  [*+let,T1  e  G-'P(2)(t,Pt3)(t)eCit,'C(T,^\ 


51  + 


+  P(t)eC(t)tD[G0,G_]C2  +  (H^2,AEVT)(t)  +  (H^,2)f)(t)  . 


We  restrict  t  to  the  interval  [ 6 , T]  where  6  is  sufficiently  large,  so  that 
(4.93)  nj,2)A  :  C([«,T]  ♦  C([S,T]> 


holds.  Then  the  following  estimate  is  fulfilled: 


{ 2  )  ~  ( 2  )  ^ 

,HT  Al  (6  T1  ■  max  ,HT  Au‘  (6  Tl  *  max 
(4.94)  lulf,  .,<1  Hul,.  ,,<1 

[6,T]  [6,T] 


(1  h'2)AuI rJt  + 
T  +  [0,T] 


+  I  h‘2)Au1,.  +  I  H(2)Au«..  ,)  < 

T  0  [6,T]  T  -  [ 6 ,T] 


<  const. ( lAl jg^Tj  ♦  6'£(4n6)3|a.[S>T])  <  ^ 


for  all  T  >  6  >  6  sufficiently  large.  This  follows  from  the  estimates  for  h'1*,  h*1* 

T  +  T  - 

given  in  Chapter  2  and  from  the  estimates  for  HQ  given  in  Markowich  (1980b).  Therefore 

(I  -  HtAE)  is  invertible  on  C([6,T)  and  we  get: 

(4.95)  VT(t>  =  ^(tK,  ♦  ^°J2  *  T*(f)(t).  tf  [6,T] 

0  0 

where  T^+»  fulfill  the  equations: 


,(2)~ 


*  P,.,!') 

rO+C.T)  9  "| 

|  G  , P  (-)P  (»)e2(" 1~2<T,f-)  G 

(1) 

. 

L  e  ej 

+  (2)  (3)  ^T '  0 

(4.97)  T*°  -  h^2)REt*°  -  ♦  [G(),GJ 
and  T1i(f)  fulfills 

(4.98 )  0(f)  -  ni2>A  0(f)  -  «'2,f  . 

T  T  T  T 

By  evaluating  at  the  boundaries  the  following  block  system  is  generated: 


(4.99) 

BE  0°<1> 

T  + 

q 

bet0°(1)  ! 

0 

j  <1  ] 

8  -  BE  V(f)(1) 

S(T)ET\i»+(T) 

S(T)E  t_<T)  I 

k  j 

1 

XT)  -  S(T)ET0(f ) (T) 

Obviously,  the  matrix  in  the  (1.1)  position  is  bounded  as  T  ♦  00 .  From  (4.97)  we  derive: 

(4.100)  S(T)E  0°(T)  =  S  (T  )E$ (T )  [G.  ,G  1  +  S(T)E(Hl21AE  0°MT)  . 

T  -  0  -  T  T  - 

(2  ) 

From  the  definition  of  we  conclude  that 

T 

T 

(4.101)  S(T)E(H^,2  ’aEt*°)(T)  =  S(T)EP(T)Doe2<T!,IP  /  DQ  e'C<  S 1  s*DP_  ’  (  s )  E~  1 AE  0°  (  s )  ds  + 

6 

+  S(T)E(tH^2,At0°)  (T) 
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where  DQ  is  the  projection  onto  the  direct  sum  of  invariant  subspaces  of  D  belonging  to 

solutions  which  decay  to  zero. 

Now  from  Markowich  (1980b),  Chapter  4  we  conclude _ 

<2)~  ~  -  ( a+ 1 )  r-  — 

(4.102)  «(TH_  Au)(T)l  -  «<H_Au)(T)l  =0!T  ),  uf  C(  L 6 . T] )  . 

This  and  (B2)  guarantee  that 

(4.103)  1S(T)Et<i°(T)«  =  o(T-(a+1)(r-1  *) 
holds*  Similarly  (4.96)  implies 

(4.104)  S(T)Et+“(T)  =  S(T)EPn)  (T)  [G+,P(2)  (T)P(3)  (T)G0)  +  o(T" (a+1 1 ( r_1 ' )  . 

Markowich  (1980b)  has  proven  that 

(4.105)  -  HAEipj:  «  <HGq,GJ 
-  (2 ) 

where  H  is  defined  as  only  the  integrals  are  stretched  to  00  instead  of  T  and 

C  is  a  regular  (rQ  +  r  )  *  (rQ  +  r  )  matrix*  We  subtract  (4.97)  from  (4.105)  getting 

(4.106)  (i|l°C  -  Ti|i°)  -  HAEl)l°C  +  =  0 

or 

(4.107)  (<|i°C  -  T*°)  =  H^2)AE(<i°C  -  *°)  +  (HAEiC°C  -  H^2  ’ AEi|)°C )  ,  6  <  t  <  T 

and  therefore 

(4.108)  -A  -  /.(6/T)  <  const .  I  HA  -  H^2)Al[6(T]  . 

By  continuation  to  [1,T]  we  get 

(4.109)  lira  I -  0 

so  that 

(4.110)  *_<1)  =  <|>°(1)C  +  o(1)  as  T  ♦  -  . 

Also  we  derive 

(4.111)  S(T)ETl>(f)(T)  -  S(T)E(H^,2)f)(T)  +  o(T-(a+1)<r'1,)llfll[(5  T]  . 

Theorem  4.2  follows  now  from  (4.103),  (4.104),  (4.110),  (4.111)  by  considering  the  system 
(4.99)  as  in  the  proof  of  Theorem  4.1.  The  convergence  results  change  correspondingly  to 
(4.83),  (4.84): 

(4.112)  ly-x  I  <  const.  (T(a+1)  ( r_1 ’( »S(T)E^(T)  [G,G  ]l  +  MHf)(T)l)  +  T_E(tnT)^)  . 

T  n»T|  0  “ 
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"T 


5.  Nonlinear  Problems 


Now  we  deal  with  problems  of  the  following  form: 


y'  “  t  f(t,y),  1  <  t  <  a  e  N  , 


(5.2)  y  c  C< [1,-] )  , 

(5.3)  b( y( 1 ) )  -  0 

f  :  Rn+1  ♦  Rn  is  supposed  to  be  continuous  in  (“,y(“)).  (5.1)  and  (5.2)  imply  that 


♦•5.4)  f(“,y<*))  -  0 

holds.  So  y(“)  can  be  calculated  a  priori  as  solution  of  a  system  of  n  nonliner 


equations.  If 


. , 3f(«,y(»)\ 
rank( - ^ - J  -  n 


then  the  solution  manifold  is  discrete  so  that  the  possible  values  of  y(")  are  known  a 
priori.  This  case  has  been  treated  by  de  Hoog  and  Weiss  (1980b).  We  will  asume  that  the 


rank  of  this  matrix  is  smaller  than  n,  so  that  we  have  to  expect  a  continuous  solution 

", 

manifold  y^fp)  with  p  c  S  C  R  ,  n1  <  n.  We  assume  that  we  have  determined  such  a  n,- 
dimenaional  manifold  and  that  f ( t , • )  *  c'fR0)  for  all  t  e  [1,“]  and 


3  f  <  t ,  y  <  u  ) ) 


A(t,p) 


J  A  ( p )  t  for  t  >  t 
i=0  1 


holds.  We  calculate  the  fundamental  matrix  E(p)$(t,p)  of  the  linearized  system 

(5.7)  y'  =*  t  °A( t  ,p )  y 

as  an  asymptotic  series.  E(p)  transforms  AQ(p)  to  its  Jordan  canonical  form  JQ(p). 

Now  we  restrict  p  to  subsets  S  C  s  so  that  r+,  r^,  r_,  r  which  are  defined  for 
Jg(p)  as  in  the  last  chapters  are  independent  of  p  in  S.  Moreover  we  require  that 
there  is  a  n  *  rQ  -  projection  like  matrix  gq  independent  of  pcs  so  that 

.  - (a+1 ) r-e  (p ) 

(5.8)  l»(t,M)[G  ,G_]I  <  C(p)t  Unt)3,  e  t  ( W )  >  0 

and  that  _ 

-2 (a+ 1 ) r-e  ( p ) 

(5.9)  •f(t,ytl(p))l  <  C(p)t  ,  E2(p)  >0,  pcs 
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holds.  Assuming  that  fy  is  locally  uniformly  Lipschitz  continuous  around  y^tu). 
Markowich  (1980b)  showed  that  there  are  solutions  y  *  y(»,u,n)  in  the  space  of  functions 
in  C([4,"J)  which  decay  to  a  finite  limit  at  least  as  fast  as  t  ^a+^r  E  .  (£nt)2^ 
where  e  =  minle^.e^)  and  6  sufficiently  large.  These  solutions  fulfill  the  estimate: 

(5.10)  «  <y(t,l),h)  -  y^iw)  -  E(li)$(t,li)  [GQ,G_]nl  <  constluWfcnt)2^  (a+1)r  e 

ro+r 

for  n  *  R  .  For  many  important  applications 

(5.11)  f(t,yjy>)  =  0 

and  4> ( t , M )  [Gq,G  ]  decays  exponentially.  In  this  case  the  right  hand  side  of  (5.10) 

*  2 

contains  the  exponential  factor  ( t ,u ) [G^ ,G_] I  and  the  algebraic  and  logarithmic 

factors  change.  It  follows  from  this  analysis  that  the  boundary  value  problem  (5.1), 
(5.2),  (5.3)  is  soluble  if  the  equation 

(5.12)  .  b(y(1,u,n)J  -  0 

n  VVr- 

is  soluble  where  b  :  R  ♦  R  and  y(t,u,n)  denotes  the  continuation  to  [ 1 

(if  it  exists).  We  assume  that  b  «  C^R"). 


The  approximating  problems  have  the  form 

(5.13)  x£  -  taf(t,xT>,  1  <  t  <  T  , 

(5.14)  b(xT(1))  -  0  , 

(5.15)  .  S(xt(T),T)  -  0 

n+^  n-(n1+r(j+r_) 

and  s  :  R  +  R  .We  assume  that  we  have  obtained  a  solution 

y*  =  y(*,u*,n*)  fulfilling  (5.10)  and  that  this  solution  is  isolated,  i.e.  the  linearized 
problem 

(5.16)  w'  = 

(5.17)  w  e  c(  (1,<»)  )  , 

(5.18) 

has  only  the  trivial  solution  w  =  0. 


a  3f ( t,y*(t) ) 

w  - 


jbiXimi  w( , ,  .  o 

3y ( 1 )  '  ’ 


Using  (5.10)  we  get  for  f^  Lipschitz  continuous 


(5.19)  =  A(t,u*>  +  0(t  2  )  . 

From  Markowich  (1980b),  Chapter  4  we  derive  that  the  general  solution  of  (5.16),  (5.17)  is 
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(5.20) 


o  Vr- 

w  =  Elu*)*  (t,lJ*,0*)5.  e  e  c 


ro  =  ni  *  ro 


i|)  is  a  n  x  ( rQ  +  r  )  matrix. 

For  the  following  we  asume  that 
(5.21) 

holds,  which  means  that  the  nonlinear  problem  and  the  linearized  problem  have  solution 
manifolds  of  the  same  dimension.  The  isolatedness  of  y*  now  implies  that  the 
(r^  +  r  )  x  (rQ  +  r_)  matrix 


(5.22) 

Now  we  define: 

(5.23) 

(5.24) 


iMX_iLLL  E(u*)i|i°(1,u*,n«)  is  regular. 


(Fy)(t)  =  (t  y1  -  f(t,y),b(y(1)),S(y(T),T)  -  S(y(T),T))  , 

F  :  (C( [1,«]  )  n  C1 ( (1,-))  n  {y|t'V  f  C( [ 1 ,“) >},<•» ( ,  j )  * 

*  (C(t1,»))  x  Rn,l(v,a)l  *  lvl(1  ^  +  lal) 


and 

(5.25) 

(5.26) 


<F  x)(t) 
T 


(t  x 1  -  f(t,x),b(x(1)),S(x(T),T))  , 


Ft  :  (C  (M.TD.Ixl  -  Ixl  (1>T)  +  T  lx'l(1<T))  ♦ 
♦  (C(  [1  ,T]  )  x  Rn,l(w,b)ll 


llwll 


1 ,  T) 


Hbl  ) 


All  involved  spaces  are  linear  normed  spaces  and  the  space  on  which  F^,  is  defined  is  a 
Banach  space. ***  We  calculate  the  Frechet-derivative  of  FT ( y* )  where  y*  is  an 
(isolated)  solution  of  F(y#)  =  0  assuming  that  S  e  C^(Rn): 

<  t f  1  ( v* > ) z ;  ( t )  =  ft'V  -  z  iilxlillilh  . 

((FT(y  z  3y  Z'  3y(  1 )  '  3y(T)  > 


(5.27) 


Assuming  that 


(5.28) 


—  (t,*)  is  locally  Llpschitz  continuous  in  the  1*1,  -norm 
3y  l  >  r°° J 

around  y*  uniformly  in  t  e  [1,”)  , 


(*) 


Norms  are  always  assumed  to  be  taken  in  the  appropriate  spaces. 
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3b  3s 

r  (y(1)),  r  —  (y(T))  are  locally  Lipschitz  continuous 

dy  ( i )  dy (T ) 

around  y* { 1 )  resp.  y*(T) 


hold,  we  derive 


"  FT(y2)l  *  const,  ly^  -  y2 1 


“y*  -  V^t]  +  T  ly*'  *  yi’'[i(T]  4  COnst"  1  =  1>2 


Moreover  Theorem  4.2  assures  that  the  problem 


2.  _  ta  +  taf(t) 

3y 


3b(  y*  ( 1 ) ) 


3S( v*(T) ,T) 


z(1)  =  a 


z(T )  =  y(t) 


ro+r-  n-(r0+r-) 

is  for  all  f  e  C([1,T]),  £5  e  R  ,  Y  (T )  e  R  uniquely  soluble  if 


3S<f^l,n  fulfUls  <A2>'  (B2>'  (C2> 


E(u*M>(t)  -  E(u*)*<t,M*>  =  E(u*)P(t,P*)tD(U*)ee(t'U*) 


is  the  fundamental  matrix  of  the  problem 


3f(t,y  (U*)) 


It  also  implies  that  the  unique  solution  fulfills 

(5.38)  lxl;i  T]  +  ^  Ix’l^  T)  <  const.  T(  0+1  1  r  Unr)2  3  M  f ,  B  ,Y  (  T  )  )  I 

so  that  F'(y#)  is  invertible  and 
T 

(5.39)  l(F^,(y4))‘1l  <  const .  T*  C‘+  1  ’  *  (  inT)2  3  . 

We  have  used  that  (5.19)  holds.  From  the  nonlinear  stability-consistency  concept  m 
Spijker  (1971)  we  conclude  that 

(5.40)  Ix^,  -  y*l[(  T]  +  T~aixT  -  (y*) '  I  ,  1  T)  4  const.  T*  °+  ’  ’  C  ( i  nT  )J  ]  I  S  (  y*  (T  )  ,T  )  I 


1 


S. 


(5.41)  IS(y*(T),T)1  <  p 1T_2 1 “+1 ’ 1 ( inT )'4 3 

holds,  where  is  sufficiently  small,  x^  is  a  solution  of  FT(x)  »  0  which  is  unique 

in  a  sphere  whose  center  is  the  restriction  of  y*  to  [1,T]  and  whose  radius  is  smaller 

than  p^T  'a+1*r(inT)  2^  with  p^  sufficiently  small.  This  holds  in  the 

Ixl,.  ,  +  T  “lx-1  -norm.  From  (5.40)  we  conclude 
[1.T]  11, T]  _ 

(5.42)  I x  -  y*8 ,  .  <  const.  T (°+1 1 1 ( inT)2 j IS( y* (T) ,T) I 

T  l  •  r  Tj 

if  (5.35)  and  (5.41)  holds. 

Because  of  (5.10)  it  is  sufficient  to  require  that  (Aj),  (Bj),  (C->)  hold  for  the 
3  S  3S 

matrix  —  (y^ly'J.T)  instead  for  jjj  (y(T,u*,n*) ,T) .  Moreover  (5.41)  is  fulfilled  if 

(5.43)  S(yM(P*),T)  =  0  for  T  sufficiently  large 
and 

(5.44)  l|^  (y<0(W*),T)(y(T,U*,h*)  -  y„(li*))l  -  0(T-2 {“+1 * r-e ) ,  e  >  0  . 

In  most  cases  of  physical  interest  y*(T)  converges  exponentially  so  that  (5.44)  is 

3  S 

fulfilled  automatically.  Therefore,  if  ( A2 ) ,  (Bj),  (C2>  hold  for  ( y^ ( p * ) , T )  and  if 

(5.43)  is  fulfilled,  convergence  follows  at  isolated  solutions  and  the  order  of  convergence 
can  be  estimated  by  (5.42). 

For  the  case  when  f  is  independent  of  t  Lentini  and  Keller  (1980)  have  generalized 
the  projection  condition  and  an  example  for  the  construction  of  an  appropriate  asymptotic 
boundary  condition  in  the  other  case  will  be  presented  in  Chapter  6. 
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6.  A  Case  Study 

The  problem  we  analyse  is  a  similarity  equation  for  a  combined  formed  and  free  con¬ 
vection  flow  over  a  horizontal  plate  (see  Schneider  (1979))*  iue  governing  equations  are 


(6.1) 


y' 


x 

^3 

X 

1  „  y1 

-  2  (1  +  r)y2 


} (1  +^r)y4 


xf(x,y) ,  0  <  x  <  ®  , 


1 

0  0 

0 

'  o' 

(6.2) 

0 

1  0 

0 

y(0 )  - 

-i 

0 

0  0 

i 

i  • 

(6.3) 

y  c 

C((0. 

»])  . 

From  (6.1)  we  conclude  that 

(6.4) 

y.  - 

yjw)  » 

(U.o, 

o 

o 

u  « 

and 


(6.S) 


0  0  0  0 

0  1  0  0 

df(x.y.(y)> 

0  0  0  0 
i  k 

+  ! 

X 

0  0  1  0 
„  u 

3y 

0  0  -2  ’I 

0  0  0  -4 

0  0"|  0 

0  0  0  -  “ 

L  2  J 

L  2  J 

For  this  problem  a  -  1,  r  ■  2  hold.  We  calculate: 


A1  (U) 


(6.6) 


J„  »  E_1A„E  » 
0  0 


E  »  diag( 1,1,1,-  — )  , 
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(6.7)  J^U)  »  E  A^UlE  »  A^U)  . 

Harkowich  (1980b)  calculated  an  asymptotic  expression  for  the  fundamental  matrix  ( x , p } 
of  the  system 

(6.8)  w’  -  x(J„  +  -  J,(m))w  , 


(6.9)  *(x,u) 


W*  “ 

1 

0(x'2) 

0(x-1 

-1 

X 

0(x_1) 

0(x_1 

0 

1 

1  -  x' 

10  0  0 


P(x,p) 


2  2 
x  u  x  u 

2  ‘  ~  ‘  2*  ~  '  2* 
diag(1,x,1,x  )diag(1,1,e  ,e  ). 


0( x, M ) 


Markowich  (1980b)  showed  that  the  problem  (6.1)  has  solutions  y(*«C^«5,/U)  which  fulfills 

2 

(6.10)  <y(x,5,  ,C-/M)  -  y  (l>)<  <  const.  x2e  4  2  ,  x  suff.  large 

l  2  00 

where  the  constant  depends  linearly  on  5^ , •  These  solutions  are  in  A£  q  +  { y^ ( u ) } 
where 


(6.11)  h  .  -  { u | u( x)  -  x"  ”  U(x),  0  e  CJ[0,»)>},  e  >  0  . 

E  ,U  b 

From  (6.9)  we  conclude  that 

(6.12)  Gq  -  [1,0,0,0]T,  G0  -  (0,1,0,0]T  . 

The  simplest  choice  of  S  is  a  linear  function  so  we  set 

(6.13)  S(X)  -  (s,(X),s2(X),s3(X),s4(X)]  . 

Condition  (B2)  of  Theorem  4.1  applied  to  our  problem  gives 

(6.14)  s1 (x)  =  0(x'2)  . 

We  choose  s^fx)  r  0.  Condition  (C2)  gives 

(6.15)  (s2(x)x_1)~1  =  0(x2)  . 

Therefore  any  matrix  S(x)  of  the  form 

(6.16)  S(x)  •  (0,s2(x) ,Sj(x) ,s3(x)] 
where 

(6.17)  s,(x)  ”  const.  *  0,  s,(x)  »  0(1),  s.(x)  =  0(1) 

2  3  4 

fulfills  (A2),  (B 2),  (C2).  A  natural  choice  is  the  following  asymptotic  boundary  condition 


GQ  -  [1,0,0,0]T,  G0  -  (0,1,0,0)T  . 


s^x)  *  0 ( x  *)  . 


(s2<x)x  1 )  1  =  0(X2)  . 


s. (x)  -  const.  *  0,  s,(x)  »  0(1),  s. (x)  =  0(1) 

2  3  4 


(6.7)  J^p)  -  E  A^pJE  -  A^p)  • 

Markovich  (1980b)  calculated  an  asymptotic  expression  for  the  fundamental  matrix  $>(x,p) 
of  the  system 

(6.8)  w'  »  x(JQ  +  ^  J^M))*  , 


(6.9)  $(x,p) 


1 

0( x-2 ) 

o(x_1: 

X 

0(x-1) 

o(x_1: 

0 

1 

1  -  x' 

0 

0 

-2 

X 

P(x,p) 


2  2 
x  p  x  p 

2  -  ~  -  2*  T  "  2* 

diag(1,x,1,x  )diag(1,1,e  ,e  ). 


2(x,p) 


Markovich  (1980b)  showed  that  the  problem  (6.1)  has  solutions  y(*/5..5,»P)  which  fulfills 


2 

_  _  Ji  x 

2  4  2 

•ytx.C^Cj.P)  -  y.lpll  <  const,  x  e  ,  x  suff.  large 


(6.10)  ly(x, C.j.^.P)  -  y^tpH  <  const,  x  e  ,  x  suff.  large 

where  the  constant  depends  linearly  on  These  solutions  are  in  Q  +  { y„( U ) } 

where 


(6.11)  h  „  -  { u I u( x)  -  x-4-EU(x)#  u  e  CJ[0,“))},  e  >  o  . 

From  (6#9)  we  conclude  that 

(6.12)  Gq  -  [1,0,0,0]T,  GQ  -  [0,1,0,0]T  . 

The  simplest  choice  of  S  is  a  linear  function  so  we  set 

(6.13)  S(X)  =  (s1(X),s2(X),s3(X),s4(X))  . 

Condition  (b2)  of  Theorem  4.1  applied  to  our  problem  gives 

(6.14)  s^x)  »  0(x-2)  . 

We  choose  8^(x)  =  0.  Condition  (Cj)  gives 

(6.15)  (s2(x)x  *)  1  =  0(x2)  . 

Therefore  any  matrix  s(x)  of  the  form 

(6.16)  S(x)  »  [0,s2(x) ,s3(x) ,s3(x)) 

where 

(6.17)  a2(x)  »  const.  *  0,  s^lx)  “  0(1).  s^(x)  “  0(1) 

fulfills  (Aj),  (B2),  (C2).  A  natural  choice  is  the  following  asymptotic  boundary  condition 


Gq  -  [1,0,0,0]T,  GQ  -  (0,1,0,0]T  . 


s1  (x)  *  0 (x  *  )  • 


(s2(x)x-1)-1  =  0(x2) 


82<x)  »  const.  *  0,  s^fx)  “  0(1),  s^fx)  ■  0(1) 


(6.18) 


(0,1,0,0]vx(X)  -  0 


which  assures  convergence  of  the  order 

2  . 

(6.19)  lvx  -  yl^  Tj  <  const,  x  exp(-  j - “  x](fnx) 

where  |i*  Is  the  parameter  value  of  the  actual  solution  y( •  ,w*»£*,CJ )  of  (6.1),  (6.2), 
(6.3)  which  is  assumed  to  be  isolated.  (6.19)  holds  because  of 

(6.20)  (0,1,0,0]y-(u)  =  0  for  p  t  R 
and  because  of  (5.41). 

Numerical  calculations  can  be  found  in  Schneider  (1978). 
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